Confinement of monopole field lines in a superconductor at T ^ 
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We apply the Bogoliubov-de Gennes equations to the confinement of a monopole-antimonopole 
pair in a superconductor. This is related to the problem of a quark-antiquark pair bound by a 
confining string, consisting of a colour-electric flux tube, dual to the magnetic vortex of type-II 
superconductors. We study the confinement of the field lines due to the superconducting state 
and calculate the effective potential between the two monopoles. At short distances the potential 
is Coulombic and at large distances the potential is linear, as previously determined solving the 
Ginzburg-Landau equations. The magnetic field lines and the string tension are also studied as 
a function of the temperature T. Because we take into account the explicit fermionic degrees of 
freedom, this work may open new perspectives to the breaking of chiral symmetry or to colour 
superconductivity. 
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I. INTRODUCTION 

The confining string is an important candidate to solve 
the confinement problem of quantum chromodynamics 
(QCD). Right after QCD was proposed as the theory 
of strong interactions 0, , different ideas for confine- 
ment were developed. Nielsen and Olesen 0| suggested 
that the magnetic flux tube vortex of type-II supercon- 
ductors could be applied, after a dual transformation, 
to colour-electric flux tube strings in QCD. Soon af- 
terwards Nambu p|, 't Hooft Q, and Mandelstam 
developed the thesis that the open colour-electric string 
was able to confine quark-antiquark systems, where the 
quark-antiquark pair is dual to a Dirac ,8| monopole- 
antimonopole pair. This idea was further developed by 
Baker, Ball and Zachariasen 0,^3 and co-authors, who 
developed dual QCD, and also studied systems similar 
to a monopole-antimonopole pair in a superconductor 
in the Ginzburg-Landau framework. At the onset of 
QCD, lattice QCD was also developed by Wilson [Tl| . 
who showed analytically that the strong coupling limit 
of QCD is equivalent to a string theory. The study of 
strings at the realistic transition between the strong and 
weak coupling of Lattice QCD was further explored by 
several authors. Bali 0, 0] and Polikarkov 0, 0] 
and co-authors studied in detail the string formation be- 
tween static quark-antiquark sources in Lattice QCD. Im- 
portantly, color superconductivity was also proposed to 
exist at finite baryon density by Alford, Rajagopal and 
Wilczek [H|. 

It is clear that QCD differs from type-II supercon- 
ductors in many details. Baker, Ball and Zachariasen 
found that dual QCD is at the frontier between type-I 
and type-II superconductivity. The spin dependence of 



'Electronic address: mjdcc@cftp.ist.utl.pt 
^Electronic address: bicudo@ist.utl.pt 
t Electronic address: pdss@cfif.ist.utl.pt 



the quark-antiquark interaction, unparallel to monopole- 
antimonopole interaction, was also addressed by De Ru- 
jula, Georgi and Glashow, and 01 Henriques, Kellett 
and Moorhouse [l8| identified the role of the confin- 
ing string in the spin dependent interactions. Moreover 
Takahashi, Suganuma and co-authors [l^. |20| also stud- 
ied the three-legged string formation between three static 
quark sources, proposed to exist in Baryons. It is clear 
that this needs a SU(3) symmetry, unparallel to type- 
I or type-II su per conductors. Okiharu, Suganuma and 
Takahashi [2l], [23], further studied the tetraquark and 
pentaquark potentials. The idea of bag model, separat- 
ing the interior of the flux tube from the remaining of 
the universe, was developed by Chodos, Jaffe, Johnson, 
Thorn and Weisskopf [2jj , and this lead Jaffe and John- 
son to propose the existence of exotic hadrons, including 
glueballs, needing also a SU(3) symmetry, p-if . 

Moreover, many problems in QCD confinement with 
strings remain to be solved. For instance in the limit 
of infinitely thin relativistic strings |25| , the world sheet 
swept out by a string in space time can only be quantized 
in 26 dimensions. For example light quarks, which con- 
dense the vacuum with chiral symmetry spontaneously 
breaking 3 P0 scalar quark-antiquark pairs described by 
Bicudo and Ribeiro [2g, are not yet fully compatible 
with string confinement. These problems, and other open 
problems of QCD confinement, motivate very interesting 
and active investigations. 

Nevertheless the main inspiration of QCD confinement, 
the monopole-antimonopole pair in a superconductor, re- 
mains to be fully explored. Ball and Caticha |27|. already 
studied the monopole-antimonopole pair in Ginzburg- 
Landau type-II superconductor. The confinement of 
magnetic field lines in superconductors occurs through 
an Anderson-Higgs mechanism such that the photons of 
the electromagnetic field acquire a mass and are there- 
fore exponentially damped in the superconductor. This 
constitutes the Meissner effect discovered experimentally 
in 1933. Indeed at the boundary between the exterior 
and the superconductor, the parallel component of the 
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FIG. 1: The magnetic field when the distance d between the 
open ends of the two half strings is respectively much smaller 
and much larger than A. 



external magnetic field is damped inside the supercon- 
ductor on a scale A, the penetration length. When the 
distance d between the monopole and the antimonopolc 
is smaller than A, the Coulomb interaction is dominant. 
When d » A, the linear interaction, characterized by 
the string tension a, is preponderant. This is schemat- 
ically shown in Fig. ^ Another relevant scale in a su- 
perconductor is the coherence length, defined as the dis- 
tance over which the superconducting order parameter 
(pair density) acquires its bulk value. The ratio of these 
two lengths determines the magnetic properties of the su- 
perconductor: for k = A/£ < 1/V2 the superconductor 
is of type-I, a complete Meissner effect is observed and 
the magnetic field is completely expelled from the super- 
conductor (perfect diamagnet) over a distance defined by 
the penetration length. If k = A/£ > l/\/2 the supercon- 
ductor is of type-II and the magnetic field penetrates in 
the material through quantized vortices, whose density 
increases until their cores (with size given by the coher- 
ence length) overlap and the system becomes normal. 

In this paper we apply the Bogoliubov-de Gennes 
(BdG) equations, enabling a microscopic understanding 
of the role of electronic states in the confinement of mag- 
netic and electric fields in superconductors. This is an 
inspiring model for QCD, where colour superconductiv- 
ity and chiral symmetry breaking can both be addressed 
at the microscopic level of fermion (quark or antiquark) 
pairing. We will be interested in the situation where two 
solenoids or magnetic rods are inserted in a superconduc- 



tor. In the limit of a very thin solenoid, the open end of 
the solenoid is equivalent to a magnetic monopole plus 
a string. It is similar to a magnetic monopole and its 
associated Dirac string introduced to satisfy Maxwell's 
equation V • B = 0. This illustrates the presence of a 
heavy quark-antiquark pair in QCD. The BdG frame- 
work, applied to the monopole-antimonopole pair is de- 
tailed in Section II. In section III we show the results 
at T = 0, and compare them with the Ball-Caticha pa- 
per. In Section IV we show the results at T 7^ 0. The 
finite temperature and finite fermion density results of 
this paper are interesting to illustrate what may happen 
in heavy ion collisions or in dense stars. In Section V we 
conclude. 



II. BOGOLIUBOV DE GENNES EQUATIONS 
AND THE MONOPOLE FIELD 



A. Fermion pairing 

Consider a fermionic system with an effective attrac- 
tive interaction between the particles defined by an 
Hamiltonian 



H = |rfr^ct(r,a) (p ^ c(r,a) 

(7 

~ \ 9 I d v ^2c\Y,a)c\r,a')c{Y,a')c{Y,a) (1) 
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where c'(r, a) creates a fermion at site r with spin a. 
A is the vector potential. In BCS theory, this leads 
to fermion-fermion and antifermion-antifcrmion pairing. 
Defining the field of a pair as A(r), and considering a 
conventional s-wave superconductor, the Hamiltonian is 
decoupled like 



dr c^(r, a) 



2m 



E, 



c(v,a) (2) 



+ J rfr[A(r) C t(r,T)ct(r,|) + A*(r) C (r,l) C (r,t)] 

Here Ep is the Fermi energy. The Hamiltonian is diago- 
nalized through a canonical transformation 



c ( r 'T) = ( M «( r )^T -<( r hL 

i 



(3) 



Here the fermionic operators 7i. CT annihilate a quasipar- 
ticle in the level i and with spin a. Defining the diago- 
nalized Hamiltonian as 



H = E g +Y. E d,a^ 



(4) 
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where E g is the groundstate energy, the physical meaning 
of the new operators is that they create the excitations 
of the system, and therefore we must restrict ourselves 
to the positive energies. Using the commutation rela- 
tions of the operators with the Hamiltonian leads to the 
Bogoliubov-de Gennes equations (BdG) [2£| for the am- 
plitudes u(r) and v(r): 



'2m 



(p - -A) 2 - E F ) Ui (r) + A(r)«i(r) = E iUi (r) (5) 



-(^(p+^A) 2 -^)^(r) + A*(r) U ,(r) = E iVi (r) (6) 

A(r) is the vector potential, A(r) is the pairing function 
given by 

A(r)= fl Yl «i(rK(r)(l-2/(f? i )). (7) 

< Ei < h v d 



Here f(Ei) is the Fermi-Dirac distribution given by 

1 



E, 
k B T 



(8) 
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where T is the temperature, assumed smaller than the 
critical temperature T c , below which superconductivity 
arises. The pairing between the electrons occurs on an 
energy scale fkojj, called the Debye energy. In conven- 
tional superconductors this is the cutoff energy provided 
by the phonon exchange between electrons. An example 
of electronic wavefunctions u and v obtained in this pa- 
per is depicted in Fig. [3 The vector potential is given 
by Maxwell's equations 

47T 

VxB = VxVxA = —Jtotai (9) 

c 

which, in the Coulomb gauge ( V.A = ), is given by 



V 2 A = -^J 



total 



(10) 



The current density originated in the supercurrents is 
obtained self-consistently by 

J(r) = ^£/(£ iK (r)(V-^A(r)H(r) (11) 
im * — 4 Tic 

i 

+ (1 - ./-(^)K(r)(V - £A(r))<(r) - c.c. 

Let us take the order parameter in the form 

A(r) = A(p, z)e~ inv (12) 

This form describes a magnetic flux equal to n flux 
quanta ($ = n(f> a — The wave functions itj and 

Vi are expanded in a way similar to ref. [29| 

,i(n-n/2)ip ^3) 
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FIG. 2: The p dependence, in the equatorial plane {z = 0), 
of the lowest electronic wavefunction. Notice that this is a 
localized wavefunction. We show in a) wavefunction u, and 
in b) the wavefunction v 



Vi(r) 



dv,3,k<Pj,IJ.+n/2,ke ^ ' ' V 



where 



V2 



III j 



- J ri 



(14) 



(15) 



Here p is an half-odd integer if n is odd and an integer 
if n is even. J m is a Bessel function. The system is 
placed in a cylinder of radius R and height h. Given the 
azimuthal symmetry of the system A does not depend on 
ip; therefore the Hamiltonian of the BdG equations may 
be simultaneously diagonalized for each value of p. It is 
therefore enough to diagonalize the matrix 



T- A 



(16) 



where 



jj'kk' 



2m J 



M±n/2 



h 2 



)Sjj'5kk' 



:Il — Ei 



(17) 



with 



h = 



dz I pdp 

I! JO 



<Pj,p±n/2,k(p: z ) 



AJp,z) 



<Pj',li±n/2,k'(j>,z) (18) 
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FIG. 3: Two choices of string orientation 

fh rR 

h = / dz pdp 



<f>j,V.±n/2,k(P> Z ) A <p(P> Z ) 2( l ) 3',H±n/2,k>(P,z) (19) 

We also obtain that 

Ajj>kk> = I dz I pdp 

JQ JO 

<Pj,H-n/2,k(P> z )Mp, Z)(j> r ^ +n /2,k' (P; z) 

(20) 

B. Monopolar gauge field 

The vector potential has two contributions one due to 
the magnetic monopoles, A ext and the other due to the 
currents that are formed in the superconductor , a. Due 
to the linearity of the Maxwell equations we have that 

A = A erf +a (21) 

The contribution from the magnetic monopoles corre- 
sponding to a magnetic field 



B 



9m r 

4/7T r 3 



(22) 



where gu is the magnetic charge and r the distance to 
the monopole, is given by 



j^ext 



9m 



47r r(r — z) 



(23) 



This form has a singularity at the origin and on the line 
(0, 0, z) with z > 0. We consider two magnetic monopoles 
one located at ri = (0,0, z\), and the other located at 
r 2 = (0, 0, z 2 ). We choose z 2 > Z\, and therefore we have 
singularities at z > z 2 and z < Z\. The potential is reg- 
ularized taking r — > R = \/r 2 + S 2 . This is equivalent 



to replacing the magnetic monopole by a semi-infinite 
infinitely thin solenoid corresponding to the magnetic 
charge and the Dirac string that carries the magnetic 
flux from the infinity to the charge. The string guaran- 
tees that VB = 0. We are interested in a situation where 
the two solenoids (or magnetic whiskers) inserted in the 
superconductor originate from one monopole flux lines 
that penetrate the superconductor and are recovered in 
the other monopole. In order for the magnetic field lines 
to diverge from one monopole and converge in the other 
the magnetic field created by the two monopoles has the 
form 



9m 

47T 



P 



P 



ri[n - (z - Zx)] r 2 [r 2 + (z - z 2 )] 



(24) 

with n = yjp 2 + {z - Zi) 2 . 

Close to each Dirac string and far from the monopoles 
the external vector potential is approximately given by a 
term like 



he 

— u-ext i I 



(25) 



since we can relate gu with the total flux $ coming 
out from a monopole where $ = n ext $o. Then it can 
be shown [3(J that the BdG equations can be recast 
in a form such that they only depend on the difference 
n a — n — n ext between the vorticity of the gap function 
(n) and the vorticity of the external line (n ex t). In the 
BdG equations we just have to replace n — > n a and the 
effect of the external infinitely thin solenoid is taken into 
account exactly if n ext is an integer. In our case the 
vector potential is more complex. 

The choice of form for the vector potential implies 
that the Dirac strings go from each monopole to infin- 
ity aligned with the z-axis in such a way that the di- 
rection of flux goes from — oo to +oo. Assuming that 
gM — we have an unit of flux along the z-axis, as 
shown in Fig. [21 The flux lines diverge (converge) from 
(to) the south (north) monopoles. The flux due to the 
field lines coming out (or in) of each monopole is equally 
divided 1/2 (—1/2) between the north and the south. 
In addition we have to take into account the flux car- 
ried by the strings. Therefore considering any plane per- 
pendicular to the z-axis, the flux traversing it will con- 
tain a flux quantum. In the south part the contribu- 
tions add like: string (+1 flux), south monopole (—1/2) 
and north monopole (+1/2). In the region between the 
monopoles where a flux tube confined by the supercon- 
ductor will appear we have: south monopole (+1/2) and 
north monopole (+1/2). Notice that with this choice 
there is no string traversing an horizontal plane when 
z\ < z < z 2 . In the north part the contributions are 
such: string (+1 flux), south monopole (+1/2) and north 
monopole (—1/2). Therefore it always adds to a unit of 
flux. 

We can as well consider that the Dirac string has a 
finite extension and goes from the north monopole to the 
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south monopole, as shown as well in Fig. |3| In this 
case it is easy to see that the total flux is zero. It seems 
therefore natural to look for solutions where in the first 
case n = 1 and in the second case n — 0, reflecting the 
different total magnetic fluxes. 

In the first case it is easy to see that the effect of the 
strings is such that, in an approximate way, in the region 
z < z\ close to the string and far from the monopole, 
the BdG equations are approximately described by the 
equations where the effect of the external vector poten- 
tial is equivalent to n a = 0, since we choose n = 1 and 
the singular contribution from the string is described by 
next = 1- A similar analysis for z > z-i leads to the 
same conclusion. Between the monopoles a finite width 
flux tube forms. As shown in ref. [3(J these field lines 
are not singular and the vorticity of the equations is 
not affected: there is no string in this region. There- 
fore n a — n — n ext = 1—0 = 1. A similar analysis for the 
second string configuration leads to a situation where for 
z < z\ and z > z^ we have an effective n a = and be- 
tween the monopoles we have n a — 1. In this case there is 
a string between the monopoles (with n ext = — 1). There- 
fore the two configurations are equivalent as expected @ . 




C. Numerical method 

The considerations above only apply to the singular 
part of the vector potential due to the strings. The sin- 
gular part effectively changes the set of basis functions 
with the replacement /i ± n/2 to p ± n a /2. Therefore 
to take advantage of this basis we just add and subtract 
to the external vector potential the limiting singular ex- 
pressions due to the strings. The remaining part of the 
vector potential is regular. It has been shown by several 
authors that in this case the contribution of the terms 
involving the vector potential in the BdG equations is 
negligible or at least very small. In the following we ne- 
glect I\ and li. We may therefore divide the system into 
two regions, using the inversion symmetry around z = 0. 
For \z\ < \zi\ we solve the BdG equations taking a basis 
with n a — 1 (equivalent to the problem of a singly quan- 
tized vortex line with no external vector potential) and 
for | jar j > \zi\ we use a basis with n a — 0. We should note 
that, as shown in ref. [28], we have to consider large sys- 
tems and consequently a large basis of functions has to 
be used to properly obtain the exponential decay of the 
various physical quantities. Otherwise, the exponential 
regime is not reached. This implies that solving the full 
3d problem requires very large matrices to be diagonal- 
ized. Therefore we have to resort to this approximation 
for the z dependence. 

In this approximation the eigenfunctions of the BdG 
Hamiltonian in each region are plane waves along z. 
Therefore the potential A does not depend on z and the 
BdG Hamiltonian is independent of z. It may be block 
diagonalized in terms of the momentum along the z di- 
rection (indexed by k). In this way the components of 




FIG. 4: (colour online) Density plot for |B| for d = 480. The 
panels correspond from left to right to no superconductor in 
an external field, the total field for n = 1 and for n — 0, 
respectively. 



the Hamiltonian are reduced to 

v2 



33 



±- 



O 



"2m 



J >±Ttq/2 

R 2 



2^2 



4fc 2 7T 



SwtEf (26) 



A 33> = 



( Pj,fi-n a /2 

{p) P dp (27) 



It is important to note that the symmetry of the BdG 
equations Ui(r) — > u*(r), Vi(r) — * — u*(r) and Ei — > — Ei 
allows to reduce the solution to the positive values of [i. 
We obtain the eigenvectors and eigenvalues for negative 
values of /i using the above symmetry. 

We may find the internal vector potential solving Pois- 
son's equation V 2 a = — ^J s . Defining a v = F(p)/p, 
Posisson's equation reduces in cylindrical coordinates to 



d 2 F 



ldF 
p dp 



dp 2 p dp dz 2 
Decomposing the current as 
4tt 



d 2 F 4tt 

= J v p 

c 



J v p = K{p, z) + f3(p, z)F(p, z) 



(28) 



(29) 
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FIG. 5: Magnetic field profile for the total field (left panels) and the internal field (right panels). We consider the cases of 
n = 1 (upper panels) and n = (lower panels) . Please note that the magnitude of the field is not visible here when it is smaller 
than the size of the arrows. This information is contained in Fig. 0] The size of the arrows is chosen for better visualization of 
the field direction at each point. 



with 

K(p, z) V /(£,)N 2 (m - n a /2) (30) 

c t—*' 

- (1 - /(^))H 2 (M + + (3(p) (A ext -a)p, 

m = ~^ y Em)\u i \ a +(i-m))\v i \ a , m 

i 

we get 

d 2 F ldF d 2 F 

w--p^ + ^ =KM+mFiPh (32) 

These equations may be solved as described in ref. |30l |. 



III. RESULTS FOR T — 

The parameters used are ef = 0.5, ojd = 0.25, g/ir = 1 
in atomic units. We should note that the relevant param- 
eter is the ratio to the critical temperature. For these pa- 
rameters the critical temperature is of the order of 0.06. 
We consider the two monopoles in a system with height 
960 and radius 480. The distance between the monopoles, 
d, is varied. Let us consider that the strings carry a quan- 
tum of flux. 

In Fig. 0] we compare the absolute value of the total 
magnetic field for the cases of no superconductor, and 
in the superconducting phase considering a n = 1 solu- 
tion and comparing to the case of n — 0. The case of 
no-superconductor is the usual magnetic dipole. Notice 
that the n = 1 and the n = cases also correspond to 
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two different configurations of confinement. The mag- 
netic vortex provides a linear confining potential. In 
the n = 1 configuration the monopole is attracted by 
the antimonopole, this is the standard picture of confine- 
ment, say in a mesonic quark- antiquark system in QCD. 
In the n = configuration the monopole is attracted 
by the boundary of the superconductor, not by the an- 
timonopole. The n = case simulates a situation where 
the magnetic solenoids are being inserted from the border 
of the material. When the solenoids are far enough, it is 
favourable for the flux lines to emerge from the material. 
As the monopoles get close together it is more favourable 
to establish a n = 1 flux tube between the monopoles. 

Taking the solution of n = 1 we look for a solution 
where the flux is quantized and equal to a quantum of 
flux. Far from either monopole and close to the strings 
the magnetic field is small. In the case of n = there 
is no flux through the system and the flux lines leave 
the superconductor from the monopoles to the border of 
the system. In Fig. [3] we show the total and internal 
magnetic fields for both n = I and n = 0. This figure 
provides complementary information with respect to Fig. 

In the case of n = 1 the presence of the flux tube 
between the monopoles is evident. In the case of n = 
the field lines avoid the central region except along the 
borders, as expected. The diamagnetic behaviour of the 
internal field, obtained subtracting the external field from 
the total field is clearly displayed in the right panels of 
Fig. El In the case of n = 1 the internal field close to 
the monopoles is not strong enough to compensate the 
external field. Note however that along the flux tube 
the internal field has the same orientation as the total 
field. Interestingly the internal field in the vicinity of the 
z = region and far from the axis winds clockwise and 
counter-clockwise. In the case of n = once again along 
the flux tubes the internal field has the same orientation 
as the total field. Between the monopoles the internal 
field opposes the total field, as expected. 

Clearly it is interesting, in the analogy of the 
Anderson-Higgs mechanism to the confinement of quarks, 
to determine the confinement potential between the two 
monopoles. Taking the analogy of the monopoles with 
the quarks in QCD the fact that these are linearly con- 
fined at large distances should translate in the supercon- 
ductor to a linear potential energy when the monopoles 
are moved apart. Indeed in Fig. 0we show the magnetic 
energy of the system, 



U = 



IBI 



-dV 



(33) 



as a function of the distance between the monopoles. The 
behaviour is similar to the one proposed for the confine- 
ment problem in QCD. We show the difference in energy 
between the normal system and the superconductor as a 
function of d. 

Interestingly, the solenoid- antisolenoid magnetic en- 
ergy includes a subtle point. To illustrate it let us now 
perform an analytical computation of the energy of the 
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FIG. 6: (colour online) a) Comparison of the total magnetic 
energy in the normal phase and in the superconducting phase. 
Also, we present the energy difference between the supercon- 
ducting phase and the normal phase, b) Energy difference 
between the superconductor and the normal phase as a func- 
tion of the distance between the monopoles plus the Coulomb 
term. 



solenoid- antisolenoid pair in the normal (with no super- 
conductor) phase. The magnetic field B can be decom- 
posed in the sum of the field created by the north anti- 
solenoid 1 and the south solenoid 2. In the limit of in- 
finitely thin solenoids, each with a quantum of magnetic 
flux, the respective fields are equivalent to a monopole 
field plus a string field, 



Bi = 
B 2 = 



47r|r _ rQ |2 

+9m 

47r|r + r | 



r-ro- g m S(x)S(y)9(z - z )(-e z ) , 

-^r + r + g m S(x)8(y)0(-z - z )(e z ) . 

(34) 



The B\ + B\ contribution is infinite but it is also con- 
stant and independent of the distance. The relevant non- 
constant term comes from the mixed term. Then we get 
two different contributions for the energy. The energy of 
the mixed monopole-antimonopole term is similar to the 
usual electrostatic Coulomb energy, 



E 



U7: 



inono— antimono 



An\2z \ 



(35) 
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FIG. 7: (colour online) Effect of temperature on the magni- 
tude of the total magnetic field for n = 1. The temperatures 
are, from left to right, T = 10" 4 , 0.02, 0.03, 0.04, respectively. 
For the density scale see Fig. [I] 



FIG. 8: (colour online) Effect of temperature on the magni- 
tude of the total magnetic field for n — 0. The temperatures 
are, from left to right, T = 10" 4 , 0.02, 0.03, 0.04, respectively. 
For the density scale see Fig. [I] 



However the energy also includes the monopole-string 

terms ^lmonopole ' ^2string ~t~ ^lstring ' ^2monopole- All 

inspection of the directions of the magnetic fields shows 
that these two terms don't cancel, in fact they contribute 
equally. We get, 



2E. 



string — antimono 



dz ■ 



+9n 



47r(z + z y 



:(e- z ■ e z )g„ 



2 9m 

Air 



fin 



■1 



z + z 



4tt|2zo| ' 



(36) 



Summing the two contributions the total energy is, 



phase (no superconductor) and in the superconductor 
phase. Although our solenoids are thinner than both the 
coherence length and the penetration length, they have 
a finite width, and we are able to account for the in- 
side (string-like) magnetic field created by the solenoids. 
Notice that the - component of the magnetic energy is 
decreasing. To isolate the linear potential we also show 
in Fig. a) the difference between the magnetic energy 
in the two phases. To get the monopolc-antimonopole 
interaction, independently of the string internal energy, 
we also add to the linear interaction of Fig. a) the 
monopole-antimonopole Coulomb interaction of eq. I|35|l . 
This produces the potential depicted in b), similar to 
the linear+Coulomb potential of QCD. 



U 



4ir\2z \ 



const. 



(37) 



Notice that the final sign in cq. I|37|) is opposed to the 
normal coulomb energy of eq. (|35|l . This occurs because 
any movement of the solenoid in the magnetic field of 
the other solenoid affects its magnetization and there- 
fore changes its internal energy. In the present case, the 
magnetization of both solenoids has to be maintained 
constant, and this requires an extra energy source, say a 
battery, thus affecting the energy balance of the system. 
Nevertheless the solenoid-antisolenoid magnetic force is 
attractive, and quite similar to the charge-anticharge 
electrostatic force. 

Indeed in Fig. a) we depict the magnetic energy 
of the solenoid-antisolenoid system both in the normal 



IV. RESULTS FOR T/0 



In Figs. and 03 we study the influence of tempera- 
ture on the total field. As the temperature increases the 
flux tube between the two monopoles narrows in the in- 
termediate region and close to the critical temperature, 
where the superconducting order vanishes, it approaches 
the normal phase regime where close to the monopoles 
the field intensity extends away from the axis. As the 
temperature increases the flux tubes narrow considerably 
until it breaks down when the superconductor becomes 
unstable and the system becomes normal. This happens 
both for n = 1 and n = 0. Again, the n — case simu- 
lates a situation where the magnetic solenoids are being 
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FIG. 9: (colour online) Magnetic field, flux and current in the "tropics" (z = 360) above the north antimonopole, and in the 
equatorial plane (z = 0). The north antimonopole is at z — 240. 



inserted from the border of the material such that it is 
favourable for the flux lines to emerge from the mate- 
rial. As the monopoles get closer it is more favourable to 
establish a flux tube between the monopoles. 

We may consider different planes perpendicular to the 
string axis. The flux contained in any plane perpendicu- 
lar to the strings axis is constant. In the regions far from 
the monopoles the flux is almost entirely contained in the 
strings. This is shown in Fig. The flux profile is shown 
along a plane that intercepts the string and it is clear that 
the flux reaches its asymptotic value very close to the 
string. In the region between the monopoles the flux lines 
emerge and converge in the opposite charge monopolc. In 
Fig. El it is clear that the flux grows smoothly from the 
origin, since in this case there is a considerable value of 
the magnetic field far from the axis. Also, we present 
results for B z and J v also as a function of temperature. 

In the central equatorial region (z — 0) the magnetic 
field is maximum on the axis and then decreases rapidly, 
on the scale of the penetration length. The field is com- 
pensated by the supercurrents that also decay rapidly as 
we move away from the axis. 

In a plane such that we are already on the strings, the 
field is finite and negative at the axis but increases to 
zero far from the string. There is a cancellation between 
the field of the string and the field of the monopoles, as 
explained above. The supercurrents vanish on the string, 
then go through a maximum at a finite distance from the 
axis and then decrease far away. Note that the flux is 
already saturated on the axis, then decreases and the 



asymptotic value is recovered far from the axis. 

As the temperature increases the penetration length 
increases. As a consequence both the increase and de- 
crease of the field with distance sets on a larger scale as 
the temperature increases. This is clearly shown on the 
various quantities plotted in Fig. 

Finally, we consider in Fig. HOI the temperature depen- 
dence of the string tension. Also, we show the tempera- 
ture dependence of the bulk value of the pairing function. 
We see that as the temperature increases both the pair- 
ing function and the string tension decrease, particularly 
as the critical temperature is approached. This critical 
temperature is defined as the temperature when the pair- 
ing function vanishes. Note that Fig. 1101 shows that the 
string tension is already quite small when the pairing 
function is still finite, even though small and close to the 
critical temperature. 



V. CONCLUSION 

Motivated by the superconductor model for the con- 
finement in QCD, we extended the Ball and Caticha 
work to include both the microscopic fermionic degrees 
of freedom of the superconductor and finite temperature. 
We addressed two confining configurations, the n = 1 
configuration where the monopolc is attracted by the 
antimonopole, and the n = configuration where the 
monopole is attracted by the boundary of the supercon- 
ductor. Our framework are the Bogoliubov-de Genncs 
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equations solved in the external field lines of a monopole- 
antimonopole pair, including the associated semi-infinite 
magnetic strings. 

Notice that this is not yet QCD, here we solve a con- 
densed matter system. The advantage of studying a 
solenoid- antisolenoid pair in a superconductor is the more 
detailed nature of our results. For instance, we identify 
the well known localized fermionic states (unatainable 
via the Ginzburg-Landau approach HH), and we analyze 
the dependence on the temperature of the string tension 
a, of the fermion pairing order parameter A and of the 
magnetic toplogical flux <£>. 

Nevertheless, this study is inspiring for the under- 
standing of confinement and of temperature and den sity 
effects in QCD El S iHl El E3, El El 
E3 or in Abelian QED QEll- These effects are usually 
explored with Lattice QCD, Schwinger-Dyson equations, 
or the Ginzburg-Landau equations. We note that the 
present framework can also be extended to model QCD. 



Importantly, this extension would be relevant to color su- 
perconductivity, where the quark density is large, and to 
the quark-gluon plasma, where the temperature is finite. 

Localized electron states exist in the region where the 
supercurrents provide a topological charge. As expected 
in the Bogoliubov-de Gennes framework, we obtain one 
fermionic localized state per angular momentum |3lj |. 

We find that the magnetic field in the vortex, say in 
the equatorial plane for n = 1, decreases with increasing 
temperature. In what concerns the penetration length 
and the topological magnetic flux, it is seen from the re- 
sults that the penetration length increases with temper- 
ature (as is well known), while the topological magnetic 
flux decreases. This decreases the string tension of the 
long distance monopole-antimonopole linear attraction. 
Notice that the the string tension is already quite small 
when the pairing function is still finite, even though small 
and close to the critical temperature. 
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FIG. 10: a) Order parameter A as a function of temperat 
b) String tension a as a function of temperature. 



